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Abstract 

A chiral invariant effective Lagrangian may be used to calculate 
the three-body interactions among low-energy pions and nucleons 
in terms of known parameters. This method is illustrated by the 
calculation of the pion-nucleus scattering length. 
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Recent articles 1 ' 2 have described a systematic effective Lagrangian frame- 
work for the calculation of reactions involving arbitrary numbers of nucleons 
as well as pions of low 3-momentum. To leading order in small momenta, 
the 'potential' for such reactions is given entirely by the tree graphs in which 
only two of the pions and/or nucleons interact; further, their interaction is 
calculated using the original effective chiral Lagrangian 3 , which consists of 
terms with only the minimum numbers of derivatives or pion mass factors, 
supplemented by contact interaction terms among nucleons. The corrections 
to these two-body interactions of second order in small momenta involve not 
only one-loop graphs, but also a large number of new terms 4 in the La- 
grangian with additional derivatives, so many that not much can be learned 
about pion-nucleon or nucleon-nucleon interactions in this way. Fortunately, 
these two-body interactions can instead be taken from phenomenological 
models that incorporate experimental information on nucleon-nucleon, pion- 
nucleon, and pion-pion scattering. The only remaining contributions to the 
potential of the same order in small momenta consist of graphs in which three 
particles (or two pairs of particles) interact, their interactions given by tree 
graphs calculated from the original effective chiral Lagrangian. Thus we can 
use the three-body interactions calculated in terms of known parameters from 
the original effective chiral Lagrangians together with experimental data on 
two-body scattering to calculate all corrections to the potential of first and 
second order in small momenta. 
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This method will be illustrated here in the calculation of the amplitudes 
for pion scattering on complex nuclei. But first, a reminder of some general- 
ities. 

Consider the amplitude for a process with N n nucleons and N n pions in 
the initial state and the same numbers of nucleons and pions in the final 
state, all with 3-momenta no larger than of order m n . We wish to develop 
a perturbation theory for this amplitude, based on an expansion in powers 
of the ratio of these small momenta (and the pion mass) to some momen- 
tum scale that is characteristic of quantum chromodynamics, such as m p . 
In counting the number of powers of small momenta in any given "old fash- 
ioned" (time-ordered) diagram for this process, we must distinguish between 
energy denominators of two types. Those of the first type arise from in- 
termediate states that differ from the initial and final states in the number 
of pions and/or in the pion energies, and are therefore of the order of the 
small momenta or the pion mass. The energy denominators of the second 
type arise from intermediate states that differ from the initial and final states 
only in the nucleon momenta, and are therefore much smaller, of the order of 
the nucleon kinetic energies. A given graph is called irreducible if it contains 
only energy denominators of the first type. These are graphs for which the 
initial particle lines cannot all be disconnected from the final particle lines 
by cutting through any intermediate state containing N n nucleons and either 
all the initial or all the final pions. We shall consider disconnected as well 
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as connected irreducible graphs, because a general connected graph is built 
up from a sequence of both disconnected and connected irreducible graphs 
interleaved with small energy denominators of the second type. (However 
in all graphs considered here, each of the initial particle lines must be con- 
nected to one or more of the final particle lines, and vice versa.) This sum of 
disconnected and connected irreducible graphs is what was referred to above 
as the potential. 

Because irreducible graphs do not contain anomalously small energy de- 
nominators of the order of nucleon kinetic energies, it is easy to count the 
number v of powers of small momenta or pion masses in these graphs. For an 
irreducible graph with V$ vertices of type i, L loops, and C separate connected 
pieces, the number of powers of small momenta or pion masses is 1,2 

v = 4 - N n - 2C + 2L + Vi A, , (1) 

% 

where Aj is an index for an interaction of type i, given in terms of the number 
rii of nucleon field factors and the number di of derivatives (or powers of pion 
mass) in the interaction, by 

A t = ^ + X -m - 2 . (2) 

(In deriving this result, we count —3 powers of small momenta for each 
line passing without interaction through the diagram, because the associated 
momentum-space delta function reduces the number of momentum factors 
in the total connected amplitude by that amount.) 
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Eq. (1) is useful because chiral invariance rules out any terms in the 
Lagrangian with Aj < 0. It follows that for any given number of external 
lines, the leading irreducible graphs (those with smallest v) are the tree 
graphs (i. e., L — 0) with the maximum number C of connected parts, 
constructed solely from vertices with Aj = 0. The contribution of these 
vertices can be read off from the effective interaction Hamiltonian (in the 
interaction picture): 

#int,A=o = i( J D 2 -l)7r 2 + i( J D- 2 -l)V7r-V7r + im2( J D^-l) 7 r 2 

?-4 
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which is derived from the most general chiral-invariant Lagrangian with Aj = 
: 
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id - 2D~ 1 F^H ■ (tt x <9 tt) -m N - 2D~ l F; 1 g A t ■ (a ■ V)tt 
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— -Cs(NN)(NN) — -Ct(NctN) ■ (NaN) (4) 

where g A — 1-25 and F n ~ 190 MeV are the usual axial coupling constant 
and pion decay amplitude; t is the nucleon isospin matrix; Cs and Ct are 
constants whose values can be inferred from the singlet and triplet neutron- 
proton scattering lengths; and D = 1 + 7r 2 /F 2 . (As discussed in ref. 2, terms 



involving time-derivatives of the nucleon field are eliminated by a suitable 
redefinition of that field, while corrections to the non-relativistic treatment 
of the nucleon in (3) and (4) appear as terms in the effective Hamiltonian 
and Lagrangian with A, > 0.) The number C of connected parts is given 
its maximum value C — N n + N n — 1 by including only graphs (with one 
qualification to be discussed later) for a single irN, NN, or irir scattering, 
with all other lines passing without interaction through the diagram. 

The corrections to these leading terms with only one extra factor of small 
momenta (or m n ) arise from (a) tree graphs, with the maximum number 
C = N n +N w — 1 of connected parts, that involve a single vertex (such as those 
arising from non-zero u and d quark masses) with Aj = 1, plus any number 
of vertices with Aj = 0. The next corrections, with two extra factors of small 
momenta (or m n ), arise from (b) one-loop graphs with C = N n + N w — 1 
involving only vertices with Aj = 0; (c) tree graphs with C = N n + N w — 1 
involving either two vertices with Aj = f or one vertex with Aj = 2 (which 
serve as counterterms for the infinities encountered in one-loop graphs), as 
well as any number of vertices with Aj = ; (d) tree graphs, constructed 
entirely from vertices with Aj = 0, that have one less than the maximum 
number of connected parts, i. e., with C = N n + N n — 2. 

As already mentioned, the vertices with Aj = 2 that contribute to cor- 
rections of type (c) contain so many free parameters 4 that little of value can 
be learned by using the effective Lagrangian to calculate these corrections. 
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On the other hand, these corrections as well as the leading terms and the 
corrections of types (a) and (b) all only contribute to the maximally discon- 
nected irreducible graphs, that consist of a connected piece involving just two 
of the incoming nucleons and/or pions, plus disconnected lines passing with- 
out interaction through the diagram for all of the other incoming nucleons 
and pions. But instead of trying to use the effective Lagrangian to calculate 
such two-body interactions, we can draw on various phenomenological mod- 
els that incorporate not only chiral symmetry but the whole body of present 
experimental information about low energy nucleon-nucleon, nucleon-pion, 
and pion-pion scattering. 

There remain only the corrections of type (d), with C — N n + N n — 2. 
These are to be calculated from tree graphs involving only the Aj = Hamil- 
tonian (3), which involves no unknown parameters. These corrections arise 
from graphs that either consist of (dl) a connected piece involving just three 
of the incoming nucleons and/or pions, or (d2) two connected pieces each 
involving just two of the incoming nucleons and/or pions, plus in both cases 
disconnected lines passing without interaction through the diagram for all of 
the other incoming nucleons and pions. Graphs of type (d2) may, like the 
graphs of types (a), (b), and (c), be taken from suitable phenomenological 
models based on experimental information about two-body scattering pro- 
cesses. This leaves only the three-body graphs of type (dl), which can be 
calculated from first principles in terms of known constants. 
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Let's first see how this applies to processes involving only nucleons. Mult- 
inucleon scattering amplitudes and bound-state wave functions are found by 
solving an inhomogeneous Lippman-Schwinger or homogeneous Schrodinger 
equation with the effective potential taken as the sum of irreducible graphs. 
The graphs for the three-nucleon terms in the effective potential are shown in 
Figures (1) and (2). A cancellation (to leading order in small momenta) was 
noted in reference 2 among the graphs of Figure (1), the only graphs that 
involve the part of the effective Hamiltonian (3) that is non-linear in the 
pion field. It is instructive to look at the reason for this cancellation. These 
graphs all involve a single quadratic interaction 2F~ 1 gA_N t ■ (ir x 7r)N plus 
linear interactions of the two pion fields in this interaction with the other 
two nucleons. In each individual time-ordered graph, the time derivative in 
the quadratic interaction makes a contribution of the order of a pion energy. 
However, by summing up the old-fashioned graphs for all the time-orderings 
of these three vertices, we obtain a Feynman diagram in which energy is 
conserved at each vertex, so that the time-derivative yields a difference of 
nucleon kinetic energies, smaller by a factor at most of order m n /m N . 

This leaves the 3-nucleon graphs of Figure (2). These are genuine contri- 
butions to what we have defined as the 3-nucleon potential, but they involve 
only the contact and pion-exchange nucleon-nucleon interactions, and their 
effect is actually cancelled by terms in the expansion of the reducible three- 
nucleon graphs of Figure (3) in powers of the ratio of nucleon to pion kinetic 
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energies. Again, the reason for this cancellation is not hard to find. Al- 
though in Figure (2) we are not summing over all time orderings, so that 
these graphs do not make up a complete Feynman diagram, the sum of all 
the time-ordered graphs of Figures (2) and (3) makes up several complete 
Feynman diagrams, in which energy denominators are replaced with pion 
and nucleon propagators, and energy is conserved at each vertex. Since the 
virtual pion energies in these Feynman diagrams are equal to differences of 
nucleon kinetic energies, and hence negligible compared with the virtual pion 
3-momenta q, the pion propagators in these diagrams are just (<f + m^) _1 . 
But these Feynman diagrams with such pion propagators are just what we 
would get from the old-fashioned diagrams of Figure (3) if we were to neglect 
nucleon kinetic energies in energy denominators for states containing a pion. 
Thus we may calculate the multi-nucleon potential to second order in small 
momenta by ignoring nucleon kinetic energies in the energy denominators of 
the leading pion-exchange contributions to the potential, and ignoring the 
three (or more) - nucleon contributions altogether. This is more or less what 
nuclear physicists have always done anywayf] 

The three-body forces are more interesting in processes involving a pion. 
For definiteness, consider the low-energy elastic scattering of a pion from a 

nucleus of nucleon number A. General considerations of scattering theory tell 

*I am grateful to J. Friar for pointing out that in some treatments of the nuclear three- 
body problem the pion exchange forces are calculated neglecting nucleon kinetic energies 
in energy denominators, and that the corrections to this approximation are of the same 
order as the other corrections considered in this work. 
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us that the S-matrix element for this process is simply given by the matrix 
element between nuclear wave functions of the sum of all irreducible graphs 
with N n = A and N n = 1. (In applying the effective chiral Lagrangian 
to such processes we are making use of the fact that typical 3-momenta of 
nucleons in nuclei are of order or less.) The leading irreducible graphs 
are those in which the pion scatters off a single nucleon, evaluated using 
the Aj = vertices in the tree approximation.^ To second order in small 
momenta, the corrections to these leading terms arise from corrections to 
the pion-nucleon scattering amplitude (from loop graphs and from vertices 
with Aj = 1, 2) which can be taken from phenomonological models of pion- 
nucleon scattering, together with connected three-body interactions among 
two nucleons and the pion, calculated from tree graphs evaluated with the 
Aj = vertices in Eq. (3). The graphs for these three-body interactions are 
shown in Figure 4. 

This is a lot to calculate, but the problem becomes much simpler if we 
restrict our attention to the pion-nucleus scattering length, for which the 
incoming and outgoing pion have vanishing 3-momenta. The leading terms 
as well as the corrections to pion-nucleon scattering give a scattering length 

that (apart from reduced-mass corrections) is just the sum of the scattering 

**There are also nominally leading terms in which the incoming pion is absorbed by one 
nucleon and the outgoing pion is emitted by another, but when these are summed over 
different time-orderings they cancel. Again, this is because summing over time-orderings 
yields a Feynman diagram in which energy is conserved, but energy cannot be conserved 
in the emission or absorption of a single real pion by a single nucleon. 
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lengths on the individual nucleons. This leaves only the three-body irre- 
ducible graphs, of which the only ones that survive in the limit of vanishing 
external pion 3-momenta are those shown in Figures 4(a) to 4(f). 

It is easiest to calculate the contributions of Figures 4(a)-4(c) and 4(f) by 
noting that the sum over time orderings in graphs of each type [and lumping 
together graph 4(f), produced by the interaction term 2F~ 4 (N(t x n)N) 2 
in the Hamiltonian (3), with the other graphs] must give the same result 
as the complete Feynman diagrams of type 4(a) - 4(c) calculated from the 
Lagrangian (4) [which does not contain the interaction 2F~ 4 (N(t x ir)N) 2 .} 
The other graphs, 4(d) and 4(e), are not summed over all time-orderings 
(because the sum would include reducible as well as irreducible graphs) and 
so their contributions must be calculated using old-fashioned perturbation 
theory. These contributions to the pion-nucleon scattering length are: 
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where subscripts a, b are pion isovector indices; r, s label individual nucleons; 
q rs is the momentum transferred between nucleons r and s in their interaction 
with the pion; and t^ r ' are the Pauli spin vector and isospin vector of 
nucleon r; and (t^) ab = —ie abc is the pion isospin vector. Note that as 
a result of a partial cancellation between (6) and (7), the integrand in the 
sum of these averages vanishes for q rs — > oo, which makes the result less 
sensitive to the behaviour of the nuclear wave function at small internucleon 
separation.^] To second order in small momenta, the pion-nuclear scattering 
length is 

1 +m n /m N ^ ( r ) [4(a)] [4(6)] [4(c)] [4(d,e)] (Q s 

(r) 

where a ab is the pion scattering length on the r'th nucleon. 

This all becomes much simpler in two special cases. One is double charge- 
exchange scattering, ir + + N — > tx~ + N', where the scattering lengths a^j 
as well as the corrections (6) and (8) vanish. The other, on which we shall 
concentrate here, is pion scattering on an isoscalar nucleus. Here t^ty + 
t^Hfr may be replaced with ^5 a bt^ ■ t^ s \ and Eq. (8) vanishes. More 
important, the contributions of the nominally leading terms in the pion- 
nucleon scattering lengths vanish, because they involve an expectation value 

of J2 r ' which vanishes in any isoscalar nucleus. The first term in (9) 

^This cancellation was noted by Robilotta and Wilkin 5 in the case of pion-deuteron 
scattering. They used a different definition of the pion field, so their results for diagrams 
4(b) and 4(c) were different from (6) and (7), but the sum of their results agrees with 
what would be found for this process from the sum of (6) and (7). 
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arises only from "a-term" corrections to the pion-nucleon scattering lengths, 
and is therefore relatively small, making it feasible to compare calculations 
of the corrections considered here with experimental measurements of the 
pion-nuclear scattering lengths. 

This may be illustrated in the paradigmatic case of pion-deuteron scat- 
tering. To evaluate the two-body terms here we need to use isotopic spin 
invariance to derive the pion-neutron scattering lengths from measured val- 
ues of the ir + p and ir~p scattering lengths. This is not entirely straightfor- 
ward, because we are interested here in the relatively small corrections to 
the leading soft-pion results for which + a nn = 0, and these corrections 
arise in part from "sigma terms" proportional to u and d quark masses that 
do not even approximately conserve isospin. Fortunately to first order in 
quark masses the isospin violation in the sigma terms affects only processes 
involving at least one neutral pion, 6 so that isospin relations can be used to 
calculate a nn . This gives the two-body terms in the 7r — d scattering length 
as 7 \X^ n % N d Kp + a ^n] = -(0.021±0.006)m- 1 . Shifting to coordinate space, 
the remaining corrections are given by: 




(10) 
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where u and w are the s-wave and d-wave parts of the deuteron wave function, 
normalized so that 



The rescattering term (10) [but not (11)] has been previously considered 
in the books of Eisenberg and Koltun and Ericson and Weise. 7 Because 
of the anomalously large radius of the deuteron, this term is considerably 
larger than the remaining three-body term (11), so it should be calculated 
including first-order corrections to the pion-nucleon scattering vertices in 
Figure 4(a). Fortunately these corrections can be taken from the measured 
values of the scattering lengths. 8 In this way one finds that 7 a' 4 ^ = 
— (0.026 ± O.OOl)/?!^ 1 . The remaining three-body terms (11) are calculated^ 
to be a [4(fc ' c)] = -0.0005m- 1 (mostly arising: from the interference between s- 
wave and d-wave parts of the wave function), in agreement with the numerical 
result quoted in reference 5. This is small compared with the uncertainties in 
other terms, and so may be neglected here, though this may not be the case 
for pion scattering on heavier nuclei. This justifies the final theoretical result 
of reference 7, = —(0.050 ± O.OOQ)!^^ 1 , which is in good agreement with 
the experimental value — (O.OSGiO.OO^m" 1 . Although the use of chiral effec- 
tive Lagrangians has turned out here only to confirm previous calculations of 
pion-deuteron scattering as well as nuclear binding, the systematic counting 

of momentum factors in chiral perturbation theory has proved its value in 

$The calculation of the integrals in Eqs. (10) and (11) was carried out by R. C. 
Mastrolco and U. van Kolck, using the Bonn wave function for the deuteron. 




(12) 
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explaining (as previous calculations did not explain) just why it is correct to 
consider only certain graphs and certain terms in the effective Lagrangian. 

I am grateful for discussions with C. Dove, J. Friar, A. Gleeson, C. 
Ordonez, U. van Kolck, and J. D. Walecka. 
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